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Abstract 

o\ 

G\ We investigate the electron tunneling into the edge ol a clean weakly in- 

teracting two-dimensional electron gas. It is shown that the corresponding 
differential conductance G(V) has a cusp at zero bias, and is characterized 
by a universal slope \dG/dV\ at V = 0. This singularity originates from the 

£^ electron scattering on the Friedel oscillation caused by the boundary of the 

system. 
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It is well known that the electron-electron interaction in one dimension leads to a power- 
law singularity in the tunneling conductance at zero bias (see, e.g., [1]). When the interaction 
is weak, this anomaly may be attributed [2] to a singular backscattering of the incident 
electrons on the Hartree-Fock potential associated with the Friedel oscillation formed near 
the tunnel barrier. The Friedel oscillation affects the electrons on the Fermi surface almost 
like a periodic potential affects electrons with the wave vector on the boundary of the 
Brillouin zone. The tunneling density of states does not acquire a gap of finite width because 
unlike a strictly periodic potential, the oscillation decays with the distance from the barrier. 
However, this decay in ID is sufficiently slow (inversely proportional to the distance) to lead 
to a power-law singularity in the density of states at the Fermi level. It has been shown [3] 
that this power-law singularity survives also in a multichannel case - for tunneling into the 
edge of a long wire of a finite width. However, the corresponding exponent decreases with an 
increase in the number of the transverse modes. The question is, whether any reminiscence 
of this zero-bias anomaly persists in the limit of an infinite number of modes, i.e., in two 
dimensions. 

We show in this paper that the differential conductance for tunneling into the edge of 
a clean, interacting two-dimensional electron gas (2DEG) is singular at zero bias. The 
backscattering from the Friedel oscillation in this case does not renormalize the zero-bias 
conductance to zero, but still leads to a non-analytic behavior of G(V) at small applied 
voltages. The corresponding cusp in G(V) is characterized by a finite slope \dG/dV\ at 
V = 0, which is universal, i.e., does not depend on the interaction strength. 

To clarify the origin of the cusp in the differential conductance, we investigate first 
tunneling into the edge of 2DEG with a short-range electron-electron interaction. In this 
case, the lowest order perturbation theory in the interaction potential is applicable. Our 
calculation of conductance in a weakly interacting system is based on the Landauer [4] 
approach. Following Ref. [2] we will focus on the effect of electron-electron interactions on the 
transmission coefficient for tunneling through a one-dimensional barrier separating two semi- 
planes. We assume that the transmission coefficient of the barrier is small, |t | 2 <C 1, and 
that the barrier is homogeneous along the y-direction [5]. In the absence of interaction, the 
current I per unit length of the barrier may be obtained from a straightforward generalization 
of the Landauer formula: 

I=tT -*V)-/ r (e)] f f^o(^)! 2 , (1) 

II J-oo ZTT J-oo ZTT 

where /;(e — eV) and f r (e) are the Fermi distribution functions in the left and right semi- 
planes respectively, and e is the kinetic energy of an electron. At small bias, Eq. (I) yields 

I(V) = (2) 

where kp is the Fermi wave vector in the 2DEG, and T = {\to\ 2 ) is the bare transmission 
coefficient at the Fermi surface averaged over the directions of momenta of the incoming 
electrons. The bare transmission coefficient may depend on energy, which gives rise to the 
well-known field effect in tunneling. Therefore, the r.h.s. of Eq. (2) can be viewed as a 
linear term in expansion of the I(V) function, the next term being proportional to V 3 . 
Unlike this one-particle field effect, the electron-electron interaction leads to a stronger and 
non- analytical correction to the current (2). 
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The existence of the barrier breaks the translation invariance in the x direction, leading 
to the Friedel oscillation of the electron density. Due to this oscillation, the electron-electron 
interaction produces an additional potential, V e s(x; k), which enhances the backscattering of 
electrons. The corresponding correction to the transmission amplitude, St(k) } can be found 
in the Born approximation. To do this, we notice that the component of the electron mo- 
mentum parallel to the barrier is conserved, and, correspondingly, the electron wave function 
has the form ^^(x^y) = e tkyy tljy i (x). In the absence of the electron-electron interaction, the 
wave function ipk(x) = <f>k x (x) is the scattering state formed by the bare barrier potential. 
The correction found in the first order in the interaction potential is 

poo A 

Sifj k (x) = / g kx (x,x')V ef [(x';k)<f) kx (x')dx', (3) 



where g kx (x } x') is the Green function of the non-interacting single-electron Hamliltonian in 
the presence of the barrier. It has the following asymptotic form at x' < x } x — > oo: 



ik x (x-x') + roe -iM*+*')| x ' > 



gk ( x ,x') = { **** r 1 ' u ~ j , ~ . ~, , 

yk x \ , ) \ m _ toetkx ( x - x >) x , < q 5 \ ) 



iK k x 



where r is the reflection amplitude. For the wave incoming from x < 0, Eq. (3) gives the 
correction of the form 

Stjj k (x) ~ -^=St(k)e tkxX , x ->■ +oo, (5) 

which defines 6t(k). The non-analytic contribution to 6t(k) is determined by the Friedel 
oscillation present in the asymptotic behavior of the effective potential at large x. This 
oscillation is caused by the reflection off the barrier of all the electrons forming the Fermi 
sea, and is characteristic for both Hartree- and exchange-type terms in V e s(x; k): 

V eB (x;k) = V H (x) + V e x(x;k). (6) 

The Hartree term is local, 

V H (x) = J U ee (r - ri)ra(ri,ri)dri, (7) 
and the exchange term is given by the following integral relation: 

Vex(z;k)<M*) = / U^r ~ riH^ri)<t>kAxiy k » (y - y ^d Yl . (8) 



Here U ee {r — r\) is the electron-electron interaction potential, and the electron density matrix 
is given by 



n(r,n) = / ^n F (qW>0<U^K 9 ^- w) , (9) 



where n^(q) = 6{kp — q) is the zero-temperature Fermi distribution function, and kp is 
the Fermi wavevector. The discontinuity in the electron momentum occupation number 
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rep(q) results in the Friedel oscillation of the density ra(r, i"i), and in the effective potential 
V e ft(x; k). The latter becomes local, 



V H (x) ~ U ee {2k F ) 8 ^ sm(2k F \x\ + J7T), (10) 
7t(2&f|:e|) 4 
gp 3 

Kx(x; k) ~ U ee (0) F . sin(2^|x| + -tt), (11) 

7r(2Ki?|x|) ' 4 

for distances x from the barrier exceeding both the Fermi wavelength Xp and the range d 
of the potential U ee (r). In derivation of (10) and (11) we assumed that the transmission 
coefficient is small, |t | 2 <C 1, and the components of k satisfy the conditions: \k x — kp\ <C kp 
and k y <^ kp. The latter two conditions allowed us to express the amplitudes of oscillation in 
(10), (11) in terms of Fourier harmonics U ee (2kp) } U ee (0) of the interaction potential U ee (r). 
Substituting now U e s(x;k.) into (3), we can find the correction 6t(k X} k y ). The oscillation 
in the potential (10), (11) leads to the non-analytical at k x = kp part of 6t(k X} k y ). The 
transmission coefficient |t(k)| 2 for electrons with energies close to the Fermi level to the first 
order in the strength of the electron-electron interaction is given by 

|t(k)| 2 = |t (k)| 2 jl - a + [U ee (0) ~ U ee (2kp)\ ^^-^QiK -kp)). (12) 

Here we presented explicitly only the non-analytical part of the dependence of the transmis- 
sion coefficient on the incident wavevector. The bare transmission coefficient, |t (k)| 2 , and 
the part of the correction a ~ m[^7 ee (0) — U ee (2kp)]/fi 2 are regular at k = kp functions. 

In order to calculate the corresponding contribution to the conductance, we have to 
substitute (12) instead of |t | 2 in Eq. (1). The calculation then amounts to integrating the 
transmission coefficient over the energies that are larger than ep but smaller than ep + eV, 
where ep is the Fermi energy and eV is the applied bias. It follows from (12) that only the 
electrons incoming with the momentum k x larger than the Fermi wavevector contribute to 
the non-analytic correction. At eV <^ ep these are electrons moving in the small range of 
incident angles almost perpendicular to the barrier, and the corresponding bare transmission 
coefficient is \to(k y = 0, ep)\ 2 . Averaging Eq. (12) over all the incident angles, we obtain the 
expression for the differential conductance G(V) = dl / dV per unit length of the barrier at 
small biases: 




m 
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G(V) = G \1 - a + 7^ [M0) " U ee (2k F )\ . (13) 

The non-analytical part in (13) provides a cusp in the differential conductance at zero bias. 
The numerical coefficient 7 \to(k y = 0, ep) | 2 /(27r) 2 T accounts for the difference between 
the average transmission coefficient T that determines Go, and the relevant for the anomaly 
transmission coefficient \to(k y = 0,ep)| 2 . The value of 7 depends on the detailed shape of 
the barrier, but typically 7 ~ 1. 

We are allowed to treat the potential energy of electron-electron interaction as a small 
perturbation only under the conditions 
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mUee(O) 
h 2 



<1, 
< 1. 



(14) 
(15) 



Hence, the "bare" interaction must be weak and short-ranged. This latter condition can 
be implemented in experiments, if a metallic gate exists very close to the 2D electron gas. 
The interaction between electrons is described by the Coulomb potential only at distances 
r smaller than the separation d between the 2D electron gas and the gate. At r ^> d the 
potential U ee (r) has a dipolar asymptotic behavior, and the condition (15) is met if d is much 
smaller than the effective Bohr radius, a# = eh 2 /me 2 , where e is the dielectric constant of 
the semiconductor. Eq. (13) leads to the following estimate of the strength of the cusp: 



(16) 



In the most interesting case of a GaAs heterostructure, the thickness of a spacer separating 
the 2D electron gas from the gate normally exceeds ag. We will show that at d ^> as the 
proportional to the interaction strength factor d/as in the estimate (16) is replaced by an 
interaction-independent universal constant. 

In the limit d ^> a#, for a pure long-range Coulomb interaction, the requirement (14) 
is equivalent to the condition that the plasma parameter e 2 / ' efivp is small. This condition 
is satisfied in a sufficiently dense interacting electron gas, and the above approach correctly 
gives the leading order Hartree-type contribution to the differential conductance. However, 
the condition (15) does not hold at any density since U ee (k — > 0) diverges. Therefore the 
slope of the cusp in G(V) might be not small, and the exchange contribution to the Eq. (13) 
must be revised. Specifically, we have to go beyond the perturbative single-electron picture 
presented above and to consider the many-body effects leading to screening of the long-range 
part of the Coulomb interaction. 

To incorporate the screening effects, we use the standard [6] relation of the tunneling 
current I with the product of the densities of states, which is valid at small |t | 2 . This 



relation can be cast in the form 
the renormalized value: 



1), if one replaces the bare transmission coefficient |t | 2 by 



T(k 



\t (k y , e F )p 



7T 2 7i 4 

— — -pi(k y} e)p r (k y} 
4m / Kz 



e . 



Here 



k x (k. 



^(2me/h 2 



L2 



(17) 



:i8i 



is the x-component of the wavevector at given total energy e and fixed value of k y . The 
densities of electron states p r ^(k y} e) are related [7] to the corresponding single-electron 
Green functions G r m(x } x 1 ; k y ; e) of the isolated right and left semiplanes respectively, 



Pr(l)(k : 



Im 



d 2 



dxd. 



- ( G r m(x,x'; k y ; 



(19) 



x=x'=0 
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The latter relation for the densities of states accounts for the fact that tunneling occurs into 
the edge of 2DEG, and the electron wavefunctions should vanish at the edge of an isolated 
semiplane, 

G r (i)(x,x'; k y ; e) x=x , =Q = 0. (20) 

As we have shown above, the zero-bias anomaly in the differential conductance originates 
from the scattering of the incoming electrons on the Friedel oscillation. Although this 
oscillation is induced by the barrier, only its tail at large distances (~ h 2 kp /meV) from the 
barrier is responsible for the anomaly. In the many-body formulation the potential of Friedel 
oscillation is represented by a specific coordinate dependence of the electron self-energy E. 
The latter contributes to the Green function, 

SG(x, x>; k y -e)=j G (°>(s, x i; k y ; e)^, x 2 ; k y ; e)G^(x 2} x>; k y ; e)dx ld x 2} (21) 



and, therefore, to the density of states. It is essential, that the self-energy in Eq. (21) 
depends not only on the coordinate difference X\ — x 2} but also on the distance from the 
barrier, i.e., on the sum of the coordinates X\ + x 2 . Our aim is to investigate the induced by 
the barrier contribution to E , to show that it has the Friedel oscillation form, and to calculate 
the corresponding correction to the tunneling density of states. The outlined program based 
on the calculation of E, enables us to generalize Eq. (II) so that the long-range nature of 
the bare Coulomb interaction between the electrons is accounted for. 

As may be anticipated from Eq. (11), the lowest order in electron charge e exchange 
contribution to E is proportional to U ee (k — > 0). For the long-range Coulomb interaction, 
this Fourier harmonic diverges. The standard way to deal with this difficulty is to sum the 
most singular at k — > diagrams in each order of the self-energy expansion in U ee (k). In the 
translationally invariant system it is easy to identify these most singular diagrams, since for 
such a system the dielectric function of the electron gas /c(i"i — r 2 ; e) is diagonal in the plane- 
wave representation (see, e.g., Ref. [6]). It is well-known [8] that the described summation 
is equivalent to the replacement of U ee (k) in the formula for the leading-order contribution 
to the exchange part of the self energy by the effective screened interaction. The latter is 
given by U ee (k) / k(k; e), where /c(k; e) is a Fourier transform of /c(i"i — r 2 ; e). In our problem, 
however, the translational invariance is destroyed by the barrier. The dielectric function 
K(ri,r 2 ; e ) depends not only on ri — r 2 , but also on the distance from the barrier, i.e., on 
X\ + x 2 . Thus, the problem of diagonalization of the dielectric function becomes non-trivial, 
and the screened electron-electron interaction potential in the vicinity of the barrier can not 
be found in a closed form. 

Our task is simplified by the fact that we do not need to know this expression for the 
screened electron-electron interaction in the vicinity of the barrier. As we noted already, 
the zero-bias anomaly arises due to the scattering of the incident electrons by the Friedel 
oscillation far from the barrier, at the distances ~ h 2 kp /meV . Thus, we are interested in 
the asymptotic form of the self energy E(xi, x 2 ;k y ; e) at X\ } x 2 ^> a#, Xp. There the dielectric 
function loses all the information about the barrier and becomes a function of | X\ — x 2 \ 
only. Thus, while calculating the singular contribution to the electron self-energy, we may 
use the formulas for the dielectric function and, hence, for the screened electron-electron 
interaction in the translationally invariant case. In dense electron liquid (e 2 / ' ehvp <C 1) this 
screened electron-electron interaction is well described by the expression found in RPA: 
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U RPA ( Xl - x 2 , k y ; e)= f dq x e^-^ U f q l K \ ^ (22) 

where k(q X} k y ; e) is a two-dimensional dielectric function of the translation invariant electron 
gas, see, e. g., [9]. Below we concentrate on the leading order exchange-type term 

Y,(x 1 ,x 2 ; k y ; e)\ xi )X2>ag )Xp = -i J dq y de'UviPK{x 1 - x 2 ; k y - q y ; e - e')G (0) (x 1 , x 2 ; q y ; e ) (23) 

since, as we already noted, in the considered dense electron liquid (e 2 / ' efiv F <C 1) the Hartree- 
type contribution is correctly accounted for by the formula (13). 

It is easy to identify now the origin of the zero-bias anomaly. In Eq. (23) all the 
relevant for singularity information about the presence of the barrier is delivered by the 
Green function. Specifically, the oscillating part of the dependence of E on its coordinates 
originates in the oscillatory dependence of the Green function on the sum of the coordinates 
X\ + x 2 . It is easy to see, that far from the barrier G^°\xi, x 2 ; q y ; e) may be decomposed as 

G {0) (x 1 ,x 2 ;q y ;e) G + (x 1 + x 2 ; q y ; e) + G~(x 1 - x 2 ; q y ; e), 

and the singular contribution to the self-energy has the form: 

Y l sm9 (x 1 ,x 2 ; k y ; e) = —i J dq y dt'UnpA.(x 1 - x 2 ; k y — q y ;e — e')G + (x 1 + x 2 ; q y ; t'). (24) 

Straightforward calculations based on Eqs. (19), (21) and (24) lead to the following 
expression for the anomalous contribution to the tunneling density of states: 



where 



6p(k y - e) = Im J d dxe^V:?(x; k y ; e) , (25) 



2 ^~ 2ik FX 

Kf( x , K\ e ) = Urpa (K x - k F , k y ; e - e F ) \—— (26) 

V 7T (2k F x) ' 

wave vector K x = k x (k y} e) is given by (18), and Urpa{<\] e) = U ee (c\) / n((\] e). 

If the applied bias is small [eV < (e 2 / ' ehv F )e F ] } we may neglect the energy and wavevec- 
tor dependence of the screened interaction in Eqs. (25) and (26), and take Urpa{<\',£) 
Urpa(0; 0). Now we note that Urpa(®', 0) = fi 2 jra is independent of the interaction strength. 
In complete analogy with the problem of tunneling into a disordered interacting 2DEG [10], 
this fact leads to universality of the correction to the electron density of states. Substituting 
Eq. (26) into (25) we find for the correction to the density of states at energies close to e F : 

, ,, , 1 /2m\i / K 2 e \ 1 ' 2 j K 2 e \ , , 

6 ^ = tA-¥) K-^-V y^~ eF )' (27) 

Thus, the anomalous exchange-type contribution of electron-electron interaction to the 
transmission coefficient is given by 
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bT(k y ; e) oc \t (k y , e F )\ — = e - — — - e F ) e - — — - e F 

po(k y ;e) 7r \ Ira J \ Ira 



-2 



Sp(k y ;e) \t (k y ,e F )\ 2 ( % k y \ ( h k y 

— r - r — t F j 

28) 



where p (k y ;e) = 2rak x (k y ; e) / ttH is the density of states (19) for tunneling into an edge 
of non-interacting electron system. Substituting (28) in (1) we find the expression for the 
differential conductance to the first order on the screened electron-electron interaction: 




G(V) = G [l-a + ^ L , (29) 



where a ~ e 2 /ehv F is small and 7 ~ 1. We neglected in Eq. (29) the Hartree-type 
contribution [see Eq. (13)], since its ratio to the exchange-type contribution is of the order 
of llR F A{2k F ) /Ur F a{0) ~ e 2 /ehv F , and therefore small. Clearly, at a larger bias, eV ^ 
(e 2 / ehv F )e F} one can not approximate Urpa{<\] e) in (26) by Urpa(0, 0), since the screening 
is not effective at |q| > At these biases the exchange contribution to differential 

conductance becomes comparable with the Hartree contribution, the slope of the curve G(V) 
decreases, and the interaction-induced correction to the conductance eventually vanishes. 
The qualitative dependence of the tunneling differential conductance on the applied bias is 
shown in Fig. 1. 

For simplicity we considered above the case of the uniform barrier, where the parallel 
to the barrier component of the tunneling electron momentum is conserved. It is clear, 
however, that our result does not depend on this barrier uniformity. Indeed, the zero-bias 
anomaly is not sensitive to the changes in the bare transmittion coefficient, and originates 
exclusively from the singularity in the electron density of states p(k y ; e) at the Fermi surface. 
It is easy to check, that even if the barrier is not uniform, the tunneling current to leading 
order in electron-electron interaction is still given by the formula 

' = £ I £ - - /*» £ ,30) 

where the superscript i = r, / refers to the left or right semiplane respectively. Clearly, the 
singular behavior of the density of states 6p l (k y ;e) at the Fermi surface (at k y = 0) does 
not depend on details of tunneling. This behavior is adequately accounted for by Eq. (27). 
All the information about the shape of the barrier, i.e., about the expression for the bare 
transmission coefficient |to(k r ,k'; e)| 2 , is absorbed in the factor F(k y ; e). It is obvious, that 
F(k y ;e) is a regular function of energy and momentum at e = e F and k l y = 0. This directly 
leads to Eq. (29). 

In conclusion, we demonstrated that a weak interaction in the clean two-dimensional 
electron gas leads to a singular contribution to tunneling into the edge of electron system. 
This singularity is characterized by a non-analytical behavior of the corresponding density 
of states near the Fermi level, and leads to a characteristic universal cusp in the differential 
conductance at zero bias. 

The authors are grateful to I. L. Aleiner for critical reading of the manuscript and useful 
comments. This work was supported by NSF Grant DMR-9423244. 
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FIG. 1. The schematic dependence of the differential conductance 
G(V) = dl I dV on the applied bias V for tunneling into the edge 
of interacting two-dimensional electron gas. The amplitude of the 
interaction-induced correction is proportional to the interaction strength 
and small, a — e 2 jefiVp, and the correction vanishes at larger bias. 
However the cusp in G(V) at V = does not depend on the interac- 
tion strength. A smooth variation of G(V) due to the field effect on 
the shape of the tunneling barrier is not shown. 
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